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Abstract 

We use the 1+3 frame formalism to write down the evolution equa- 
tions for spherically symmetric models as a well-posed system of first 
order PDEs in two variables, suitable for numerical and qualitative 
analysis. 

1 Introduction 

We shall use the 1+3 frame formalism [TJ [2] to write down the evolution 
equations for spherically symmetric models as a well-posed system of first 
order PDEs in 2 variables. The formalism is particularly well-suited for 
studying perfect fluid spherically symmetric models [3], and especially for 
numerical and qualitative analysis, and is useful in various applications, such 
as structure formation in the spherically symmetric dust Lemaitre-Tolman- 
Bondi model. This preprint is intended as a resource paper for researchers 
working in this field. 
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2 Spherically symmetric models 

The metric isQ 

ds 2 = -N 2 dt 2 + ( ei l y 2 dx 2 + (e 2 2 )" 2 (^ 2 + sin 2 $ dip 2 ). (1) 
The Killing vector fields (KVF) are given by [I]: 

d v , cosf d$ — sin^cot-^ d v , simp d$ + cos cot ^ d v . (2) 



The frame vectors in coordinate form are: 



e = N l d u ei = e^d x , e 2 = e 2 2 ^, e 3 = e 3 3 c\,, (3) 

where e 3 3 = e 2 2 / sin??. iV, ei 1 and e 2 2 are functions of t and x|l 

This leads to the following restrictions on the kinematic variables: 

a a p = diag(-2(x + , a+, a+), u af3 = 0, ii a = (u h 0, 0), (4) 

where 

w 1 = e 1 lnJV; (5) 
on the spatial commutation functions: 

/ n 13 \ 
a a = (a u a 2 ,0), n a(i = , (6) 

V "13 J 

where^l 

ai = eilne 2 2 , a 2 = m 3 = ~-e 2 2 cot (7) 
and on the matter components: 

q a = (qi, 0, 0), 7i af 3 = diag(-27T + , tt + , tt+). (8) 

1 We use x instead of r because r is used to denote the spatial derivative of H. 

2 Note that the frame vectors e,2 and e3 tangent to the spheres are not group-invariant - 
the commutators [b2, d v ] and [e3, d v ] are zero, but not with the other two Killing vectors. 
The frame vectors eo and ei orthogonal to the spheres are group-invariant. 

3 The dependence of and on d is due to the fact that the chosen orthonormal 
frame is not group- invariant. However, this is not a concern, since the d dependence will 
be hidden. 
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The frame rotation Q a p is also zero. 

Furthermore, n 13 only appears in the equations together with e 2 n 13 in 
the form of the Gauss curvature of the spheres 

2 K := 2(e 2 - 2n 13 )n 13 , (9) 

which simplifies to 

2 K = (e 2 2 ) 2 . (10) 

Thus the dependence on -& is hidden in the equations. We will also use 2 K m 
place of e 2 2 . 

The spatial curvatures also simplify to: 

3 S aP = dmg{-2 3 S + 3 S + , 3 S + ), (11) 

with 3 R and 3 S + given by: 



3 R = 4eiai - 6a 2 + 2 2 K (12) 
3 S + = -|eiai + | 2 K (13) 

The Weyl curvature components simplify to: 

E a(3 = diag(-2£ + , E + , E + ), H af3 = 0, (14) 

with E + given by 

E + = Ha + + a 2 + + 3 S + - (15) 
To simplify notation, we will write 

2 K, iii, Qi 

as 

K, it, a. 

To summarize, the essential variables are 

iV,ei\ K, H, a+, a, //, q u p, tt+, (16) 
and the auxiliary variables are 

3 R, 3 S + , ii. (17) 
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So far, there are no evolution equations for N, p and 7r + , and they need to 
be specified by a temporal gauge (for N), and by a fluid model (for p and 



The evolution equations are nowQ 

eoei 1 = (-H + 2a + )e 1 1 (18) 

e K = -2(H + a+)K (19) 

e H = -H 2 - 2a\ + |(e x 2a)ix - |(/i + 3p) + |A (20) 

e o-+ = -3i7o-+ - |(ei + u + a)u - 3 S+ + tt+ (21) 

e a = (-H + 2a + )a- (e^u)^ + a + ) (22) 

e /x = —3H(p, + p) — (e x + 2-it — 2a)gi — 6a + 7r + (23) 

e gi = + 2a + )q l - e x p - (p, + p)u + 2(ei +u- 3a)ii + . (24) 

The constraint equations are the Gauss and Codazzi constraints, and the 
definition of a: 

= 3H 2 + fR - 3a 2 + - fi - A (25) 

= -2e 1 (H + a+) + 6aa + + qi (26) 

= (ei - 2a) K, (27) 

where the spatial curvatures are given by 

3 R = 4e t a - 6a 2 + 2K (28) 

3 S + = -|eia+ \K. (29) 



3 The matter and gauge 

There are various choices for the matter. 

3.1 Perfect fluid 

A perfect fluid is defined by [2] 

T ab = fiu a u b + p(g ab + u a u b ). (30) 
4 We include a non-negative A. 
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with p to be specified. In general, the 4- velocity vector u of the perfect fluid 
is not aligned with the vector e of a chosen temporal gauge. In spherically 
symmetric models, u is allowed to be of the form 

u = r(e + uei), T = {l-v 2 )- 1 *. (31) 

We choose a linear equation of state for the perfect fluid: 

p = (7 - (32) 

where 7 is a constant satisfying 1 < 7 < 2. Then we obtain for the tilted 
fluid: 

p = — 2 p (34) 

?i = ^ (35) 

where G± = 1 ± (7 — l)v 2 . Thus p, qi and 7r + are given in terms of p and v. 
These are then substituted into the evolution and constraint equations. 
The evolution equations for p and qi now give (in terms of p and v) 

e /i = -y^e lf i - ^fpe ± v - -^-p [(3 + v 2 )H + 2v(u - a) - 2v 2 a + ] (37) 

( 7 -l)(W) 2 [(3 7 -4)-( 7 -l)(4- 7 )«> 
7 G_/i ^ + G+GL 1 

- ^— — - [-(3 7 - 4># - 2va+ + G-ii + 2(7 - l)v 2 a] . (38) 
3.2 Scalar fields and anisotropic fluid 

The total energy-momentum tensor of a non-interacting scalar field (f> with a 
self-interaction potential V{4>) is 

T a f = ;a ;b - g ab {\<P-A'' c + (39) 
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where = <fi(t, x). In particular, exponential potentials have been the subject 
of much interest and arise naturally from theories of gravity such as scalar- 
tensor theories or string theory [5] . Spherically symmetric scalar field models 
have been studied in [6]. 

A spherically symmetric model can also admit an anisotropic fluid mat- 
ter source, in which the energy momentum tensor has energy density fi, a 
pressure pn parallel to the radial unit normal and a perpendicular pressure 
pj_. Fluids with an anisotropic pressure have been studied in the cosmolog- 
ical context for a number of reasons [TJ. An energy-momentum tensor of 
this form formally arises if the source consists of two perfect fluids with dis- 
tinct four- velocities, a heat conducting viscous fluid and a perfect fluid and a 
magnetic field; in addition, a cosmic string and a global monopole are of the 
form of an anisotropic fluid. Most importantly, perhaps, a contribution in 
the form of an anisotropic fluid arises when averaging the Einstein equation 
to obtain the averaged field equations in spherically symmetric geometries 

3.3 Temporal gauge 

The common temporal gauges used in spherically symmetric cosmological 
models are the synchronous gauge and the separable area gauge. The syn- 
chronous gauge is useful when used with a dust perfect fluid (7 = 1) because 
the dust perfect fluid has zero acceleration (u = 0). This gives the Lemaitre- 
Tolman-Bondi models. It may also be useful when used with a non-dust 
tilted perfect fluid. The un-normalized system is well-posed when the Gauss 
and Codazzi constraints are used to eliminate the spatial derivatives. H- 
normalization preserves well-posedness. 

The separable area gauge has a special case (a = 0), called the timelike 
area gauge. The un-normalized system is well-posed when eo(H + a + ) is 
used, and the Gauss constraint is solved for H — a + . (/J + a + )-normalization 
preserves well-posedness, while iJ-normalization does not. 

4 Special cases with extra Killing vectors 

Spherically symmetric models with more than 3 KVF are either spatially 
homogeneous or static. Let us discuss the spatially homogeneous cosmologi- 
cal models. Spatially homogeneous spherically symmetric models consist of 
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two disjoint sets of models: the Kantowski-Sachs models and the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) models. Static and self-similar spheri- 
cally symmetric models have been studied in [6], [91 [3] . 



4.1 The Kantowski-Sachs models 

The spatially homogeneous spherically symmetric models (that has 4 Killing 
vectors, the fourth being d x ) are the so-called Kantowski-Sachs models [I]. 
The metric ([1]) simplifies to 



ds 2 = -N(t) 2 dt 2 + {e 1 1 {t))- 2 dx 2 + (e 2 2 (t))- 2 (^ 2 + sin 2 ^ dip 2 )- (40) 



i.e., N, ei 1 and e 2 2 are now independent of x. 

The spatial derivative terms e x ( ) vanish and as a result a = = u. Since 
ii = 0, the temporal gauge is synchronous and we can set iV to any positive 
function of t. 

The Codazzi constraint restricts the source by 



p and 7r + are still unspecified. 

The evolution equations for Kantowski-Sachs models with unspecified 
source are: 



qi = 0. 



(41) 



e ei 



i 



-H + 2a + )e 1 1 
2(H + a + )K 
H 2 -2o\ - |(^ + 3j9) + 
3Ha + -\K + n + 
3H(fi + p) — 6<T + 7Y + 



(42) 
(43) 
(44) 
(45) 
(46) 



e K 
e H 




e a+ 



The remaining constraint equation is the Gauss constraint: 



= 3H 2 + K - 3a 2 - y, - A. 



(47) 



The spatial curvatures are given by 



3 R = 2K 



(48) 
(49) 



3 S+ = \K. 
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4.2 The FLRW models 



Spatially homogeneous spherically symmetric models, that are not Kantowski- 
Sachs, are the Friedmann-Lemaitre- Robertson- Walker (FLRW) models (with 
or without A). The source must be of the form of a comoving perfect fluid 
(or vacuum). 

The metric has the form 

ds 2 = -N(t) 2 dt 2 + £ 2 (t)dx 2 + £ 2 (t)f 2 (x)(d$ 2 + sin 2 $ dp 2 ), (50) 

with 

f(x) = sins, x, sinhx, (51) 

for closed, flat, and open FLRW models respectively. The frame coefficients 
are given by ei 1 = ^(t) and e 2 = £ _1 (t)/ _1 (x). Then a+ = |e ln(ei 1 /e2 2 ) 
vanishes. N = N(t) implies that u = 0; i.e., the temporal gauge is syn- 
chronous, and we can set TV to any positive function of t. The Hubble scalar 
H = e \n£(t) is also a function of t. [f| 

For the spatial curvatures, 3 S+ does vanish because fl5T|) implies e^a — 
while 3 R simplifies to 

3 R = ^, k = 1,0,-1, (52) 

for closed, flat, and open FLRW respectively. 

The evolution equation for a + and the Codazzi constraint then imply that 
7r + = = qi] i.e., the source is a comoving perfect fluid, with unspecified 
pressure p. 

The evolution equations simplify to: 



= m (53) 

e H = -H 2 - + 3p) + | A (54) 
e /i — — 3H(/i+p) (55) 



The Gauss constraint simplifies to 



3k 

= 3H 2 + — -fi-A, k = 1,0, -1. (56) 



5 We shall not list the KVFs as they are complicated in spherically symmetric coordi- 
nates and not needed here. 

6 That eia does not vanish is consistent with the frame vector ei not being group- 
invariant. 
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Note that /i and p also depend on t only, and that p is not specified yet. 

The vacuum cases are the de Sitter model (A > 0, k — 0), the model 
with A > 0, k — 1, the model with A > 0, k — -1, the Milne model (A = 0, 
k = — 1), and the Minkowski spacetime (A = 0, k — 0), which is also static. 
The model with A > 0, k = 1 is past asymptotic to the de Sitter model with 
negative H and is future asympotic to the de Sitter model with positive H. 
The model with A > 0, k = — 1 (and positive H) is past asymptotic to the 
Milne model and is future asympotic to the de Sitter model with positive H. 

5 Synchronous gauge, tilted perfect fluid 

We shall investigate perfect fluid models with linear equation of state us- 
ing the synchronous gauge. We shall simplify the equations step-by-step, by 
choosing the synchronous gauge, eliminating spatial derivatives, and speci- 
fying the perfect fluid. 

The equations in synchronous gauge (it = 0) are: 




(-H + 2a + )e 1 1 



(57) 
(58) 
(59) 
(60) 
(61) 
(62) 
(63) 




e a + 



e H 



e <?i 



e a 



-H 2 - 2a 2 + - + 3p) + ±A 
-3Ha + - 3 S + + 7T+ 
(-H + 2a+)a-e 1 (H + a+) 
—3H(fi + p) — (ei — 2a) qi — 6a + n + 
{-AH + 2cr + )gi - + 2{e x - 3a)vr + . 



The constraint equations are: 



= 3H 2 + - 3a 2 _ - fi- A 
= -2ei(H + <r+) + Qaa + + q 1 
= (ei - 2a) K. 



(64) 
(65) 
(66) 



where the spatial curvatures are given by 



3 R = 4e ia - 6a 2 + 2K 
*S + = -| ei a + \K. 



(67) 
(68) 
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The evolution equations (I60p and (15Tj) contain spatial derivative terms, 
but these can be replaced using the constraints (154"j) and (1551) : 

e 1 a = -|i/ 2 + |a 2 -|ir+|a2 +1^ + ^ (69) 
e 1 (F + a + ) = 3aa + + (70) 

As a result, equations (1601) and (16T1) now read: 

e a+ = -3#ct + - |# 2 + fa 2 -\K+ \o\ + ±/x + ±A + tt + (71) 
e a = -(if + o+)a- (72) 

The benefit here is that the evolution equations for the geometric part are 
now free of spatial derivative terms. 
The spatial curvatures are given by 

3 R = -6H 2 + 6or 2 + 2/i + 2A (73) 
3 S + = \H 2 - \a 2 + \K- \a\ -\/i- §A. (74) 

Lastly, we specify the perfect fluid with linear equation of state. From 
equations (132]) - (138]) and the above equations, the final form of the system is: 

eoex 1 = (-H + 2a+)e 1 1 (75) 
e K = -2(H + a + )K (76) 

e^-^-^- ^-^-^V tA (77) 



6G H 

e a + = -3Ha + - \H 2 + \a 2 - \K + \a 2 + + (1 ~ ( J + + ±A (78) 



e a = -(H + a + )a- -—/i (79) 



2G 

7 .. r^o , „.2 



/1 [(3 + v 2 )H -2va-2v 2 a + ] (80) 

eol , + fr-V-^e,, - K»7-4)-(7-l)(4-7W V 
■y&^/i G + G_ 

(1 - w 2 ) 

[-(37-4)fl"-2o- + +2(7-l)«a]v, (81) 
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where G± = 1 ± (7 — l)v 2 . The constraints are: 



eict 



\E 2 + fa 2 - \K + fa 2 



e\(H + + ) = 3aa + + 



2G_ 



-/' 



eiK = 2aK. 
The spatial curvatures are given by 



3 R 
3 S+ 



—QH 2 + 6a 2 + 2/x + 2A 

1 rr2 1 2 1 I17 12 



i/ 1 



(82) 
(83) 

(84) 



(85) 
(86) 



5.1 Well-posedness 

We now show that the system is well-posed for 7 > 1. 

The coefficient matrix for the spatial derivative terms is: 



/ 



V 



7t> 

(7-l)(l 



,2\2 



7g- 



)2 [(3 7 -4)-( 7 -l)(4- 7 )t;> 



7G-// 



G4-G- 



Its eigenvalues are 



(2- 7 )t;± Vt^T(1 



with corresponding eigenvectors (for example) 



(87) 



(88) 



/i7(l + V7 - lw) 2 



Vt^T(i - ^ 2 )g + 

//7(1 - V7 - l^) 2 
v 7 ^ 



(89) 



The matrix is diagonalizable for 7 > 1, with c s = ^7 — 1 being the speed 
of sound in the perfect fluid. The system fl75l) - fl8~T]) is thus well-posed for 
7 > 1. For 7 < 1 the system is elliptic and not well-posed. 



7 Strictly speaking, we should also include the factor ei 1 /N in the matrix, but the result 
on well-posedness is the same. 
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6 Irrotational dust (Lemaitre-Tolman-Bondi 
model) 

The Lemaitre-Tolman-Bondi (LTB) model [T01I7] is the spherically symmetric 
dust solution of the Einstein equations which can be regarded as a generaliza- 
tion of the FLRW universe. LTB metrics with dust source and a comoving 
and geodesic 4-velocity constitute a well known class of exact solutions of 
Einstein's field equations [U [7] . 

For the dust case 7 = 1 with zero vorticity, we can use the freedom 
within the synchronous gauge to set v — 0, so that the synchronous frame is 
comoving with the perfect fluid and we obtain: 



eoe/ = (-H + 2a+)e 1 1 (90) 

e K = -2(H + a + )K (91) 

e H = -H 2 - 2a\ - §// + ±A (92) 

e a + = -3Ha+ - \H 2 + \a 2 - \K + \a% + |/x + |A (93) 

e a = —(H + a + )a (94) 

eo/W = —SHfjL. (95) 



Notice that the system is completely free of spatial derivatives, and is thus 
well-posed. 

The constraints are: 

e x a = -\E 2 + fa 2 -\K+ \a 2 + + |/i + |A 
ei(H + cr + ) = 3aa + 
e t K = 2aK. 

The spatial curvatures are given by 

3 R = -6H 2 + 6a 2 + + 2/J + 2A (99) 
3 S + = \H 2 - \a 2 + \K- \o\ - i/x - \k. (100) 

Further simplifications with a and K are possible. 

A suitable normalization factor is (3 = H + <j + , introduced for G2 models 
in [21 H], and used for the LTB model in [llj . With this normalization, it can 
be shown that at late times the LTB solutions that are ever-expanding will 
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(96) 
(97) 
(98) 



tend to the isotropic and homogeneous Milne solution, with the following 
rates: 



(3 ~ e~\ 



a 2 -K 



1 ~ e 



— T 



(101) 




— T 



~ e 



— T 



(102) 



That is, the rates are the same for all dust observers, although the multi- 
plicative "constants" depend on the radius. This dependency reveals itself 
in the leading order of ratios of variables such as the density contrast. 

6.1 Structure formation 

Structure formation in the LTB model has been studied in [12]. More re- 
cently, the LTB inhomogeneous dust solutions have been examined numeri- 
cally and qualitatively as a 3-dimensional dynamical system, in terms of an 
average density parameter, (Q) (which behaves dynamically like the usual Q 
in FLRW dust spacetimes), and a shear parameter and a density contrast 
function which convey the effects of inhomogeneities |llj . The evolution 
equations for the averaged variables are formally identical to those of an 
equivalent FLRW cosmology, and are an alternative set of evolution equa- 
tions to those presented above. In particular, the phase space evolution of 
structure formation scenario was examined in [TT] . 
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